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Abstract 

r-| I This paper represents an in-depth treatment of the chain inflation scenario. We fully 

determine the evolution of the universe in the model, the necessary conditions in order 
cn ■ to have a successful inflationary period, and the matching with the observational 

results regarding the cosmological perturbations. We study in great detail, and in 



> 



Tjj- I general, the dynamics of the background, as well as the mechanism of generation of 

the perturbations. We also find an explicit formula for the spectrum of adiabatic 
perturbations. Our results prove that chain inflation is a viable model for solving 

Q ■ the horizon, entropy and flatness problem of standard cosmology and for generating 

^ . the right amount of adiabatic cosmological perturbations. The results are radically 

different from those found in previous works on the subject. Finally, we argue that 
there is a natural way to embed chain inflation into flux compactified string theory. 

/\^ • We discuss the details of the implementation and how to fit observations. 
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1 Introduction 

String theory is our best developed theory of quantum gravity and should describe 
the evolution and birth of our universe. In particular, it should provide a mechanism 
for inflation. Unfortunately, all models constructed so far suffer from various defl- 
ciencies. A signiflcant difliculty is that potentials appropriate for slow-roll seem to be 
uncommon in string theory. On the other hand, the theory is in many cases poorly 
known, and therefore so are the corrections to the potentials, which can be important 
for large values of the flelds. 

Most successful models of inflation are based on second order phase transitions 
with a fleld that rolls slowly through fleld conflguration space. In the simplest models 
the field travels long distances of the order of the Planck mass, which makes a low 
energy approximation questionable. There is also the issue of fine-tuning of the models 
if one wants to obtain enough infiation. 

In this paper we will explore the possibility of using a series of first order phase 
transitions to generate infiation. While it is well known that a single first order phase 
transition, as in the case of old infiation, cannot do the job, the situation could be 
dramatically different with a large number of transitions. The idea is called chain 
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inflation [Il[2]. The basic idea is that inflation need not proceed in one step. Instead 
there is a long series of metastable minima where the inflaton field or fields spend 
just a short time allowing for a fraction of an e-folding. The field then tunnels to the 
next metastable minima in the series. Only after adding up several such steps will 
there be a sufficient number of e-foldings. 

It is widely believed that string theory has a rich structure of many metastable 
vacua separated by potential barriers related to domain walls (D- and M-branes). 
Such a setup favours first order transitions, and should provide suitable conditions 
for chain inflation. The basic idea is that the early universe searched its way through 
the landscape experiencing a large number of tunneling events before settling down in 
the present state. During this process the universe should have gone through periods 
of inflation. 

Our motivation for studying chain inflation is twofold. On the one hand, this 
scenario has not been studied in sufficient detail to make it robust. Moreover, results 
in previous works are not fully consistent with a complete treatment of the subject, 
[D [21 El m El [6l [7] . On the other hand, our interest is directly related to a previous 
study of the topography of the string landscape. In P, E] it was discovered that series 
of simply related vacua are a generic feature of the string landscape in the context of 
flux compactifications (for a review see for example fiUi [llj and references therein). 
Contrary to previous studies of chain inflation, this provides us with a reasonably 
constrained framework for implementing chain inflation. Furthermore, in this case 
we are able to propose an inflationary model where the agents are complex structure 
moduli, in contrast to the models so far described in the literature, which involve 
Kahler moduli or open string moduli. 

The paper is divided into four main sections. In section [2] we discuss the basic 
problems of inflation with first order phase transitions, and discuss the various con- 
straints that need to be met for successfully solving the horizon, entropy and flatness 
problem and achieving percolation and large scale thermalization. In section El we 
discuss the generation of cosmological perturbations in chain inflation, proposing a 
specific set of equations for determining them and computing the power spectrum. 

In the first sections our discussion is at a general level, while a detailed model 
based on string theory is specified in section |H There, we review the construction 
of a series of minima using flux compactification and apply it to chain inflation. We 
compare our results to the experimental data in section El We then end with some 
conclusions, speculations and outlook. 



2 Successful chain inflation 



The key-point of chain inflation is that the slow-rolling of the inflaton is substituted 
by a series of subsequent tunneling processes. Chain inflation must still fulfill the 
same demands as ordinary inflation with a slow evolution of the energy density of the 



vacuum allowing for an accelerated expansion of the universe. At the same time, in 
order to obtain a homogeneous and isotropic universe, we need the tunneling events 
to proceed in such a way as to guarantee percolation and homogenization at every 
step of the ladder of transitions. 

While these requirements are enouglo to solve the horizon, entropy and flatness 
problems of standard cosmology (we can call it background requirements) , the model 
must also provide the right amount of density perturbations and the correct features 
of the spectrum. 

We will deal with this last aspect later in the paper (section [3]), for the moment let 
us proceed by establishing the background requirements for the tunneling processes. 

2.1 Basic setup 

Consider a theory with a potential with many metastable minima at different ener- 
gieqj. The number of fields could in principle be very large, as it is supposed to be 
for string moduli, and their dynamics as they pass through the various minima could 
be complex, with rolling, tunneling, and jumping phases. In the following we will 
simplify the description by modeling the dynamics purely through the stress energy 
tensors of the fields and the gravity field. 

The phase transitions among the various vacua proceed by nucleation of bubbles. 
The bubbles are filled by the new vacuum at a lower energy, and all the latent heat 
of the transition is stored within the walls of the bubbles. When bubbles collide the 
energy of the walls is converted in radiation (massless string states). The process of 
bubble collision and the transfer of the energy to radiation are complex phenomena, 
that can be studied numerically. We are not going to do so, and will instead stick to 
a simplified model based on a number of assumptions: 

• we consider values for the decay rate per unit physical volume r„ between vacua 
n and n — 1 such that percolation and large scale thermalization is achieved 
at every phase transition (we will investigate which values allow this in the 
following) 

• decays between two distant vacua are strongly suppressed 

• the bubbles are produced with negligible volume and expand at the speed of 
light 

• the energy of the walls is optimally transferred to radiation in the collision and 
the thermalization time is very short 



^A viable model for inflation must also provide a suitable mechanism for reheating, but we will 
not address this question in this paper. 

^If we are allowed to consider the string landscape as governed by a very complicated effective 
action with a potential depending on many fields, this description applies to that case as well. 



In this way, we can assume that a homogeneous and isotropic space-time is ob- 
tained in a very short time in every patch and phase, and we can therefore describe 
the metric using the Friedman- Robert son formula 

ds'^ = -dt^ + a{tfd^^. (1) 

Schematically we have the following picture. As time passes the universe is divided 
into regions of different phases (new vacua) created through bubble nucleation and 
rapid percolation. The collisions among bubble walls produce an amount of radiation 
which remains trapped within every phase. Our actual Universe is dominated by 
one of these phases. From our point of view, therefore, we observe the Universe as 
filled with the radiation produced from the collisions of the walls of the bubbles of 
our phase and small patches, distributed in an homogeneous way, composed of the 
previous phases (which are contained one within the other, together with the radiation 
produced at every step). 

It is convenient to write the equations involving the energy tensors of the vacuum 
and the radiation by formally distinguishing every phase with an index m and asso- 
ciating to each one a vacuum and a radiation stress-energy tensor, distinguished by 
£ = V,r 

where m^'*", p^'™", P^'"*, for £ = {r, y}, are the proper four- velocity, energy and pres- 
sure densities. We write the energy density, pressure and stress-energy tensors for 
radiation and vacuum as 



P 



Y^p^m p^^Y^p^m Tl^^Y.^'^-^ (3) 



m 



and the total energy density, pressure and stress-energy tensor as 

P^Hp' P^Hp' T,^^Y.n^. (4) 



We furthermore write 



P^^P'^ + P", (5) 

expressing that the "vacuum energy" is actually the sum of the energy density in the 
interior of the bubbles, p^, plus the energy stored in the bubbles walls, p^ . 

The energy and pressure densities of our generalized fluidqj depend on space and 
time, but, due to the hypothesis of homogeneity and isotropy within every phase, 
their average over the volume occupied by the phases depends only on time 

(p^-'-(x,t))=p^-'-(t). (6) 



■^As well known, the hydrodynaniical fluid description do not apply for scalar fields and dark 
energy, and therefore we cannot use it for modeling the dynamics of our vacuum fields here. Nev- 
ertheless, it is possible and well-defined to speak of energy density and radiation with reference to 
the components of the energy tensor, that is why we speak of generalized fluids. 



Written in terms of p and p'', the Friedman and Chauduri equations are 

H' = ^(P^ + /) (7) 



H = -AnG J2 (/ + P') 



£=V,r 



2.2 A slow roll through the landscape 

In order to have a successful chain inflation mechanism, the system must undergo a 
series of successive tunnelings such that the decay of the vacuum energy is effectively 
slowed down. In this way inflation can effectively cure the problems related to entropy, 
flatness and horizon size, by lasting for a sufficient number of e- foldings (A^e ~ 60). 

To avoid a too rapid decay of the vacuum energy it is necessary that each tunneling 
step occurs only after a suitable and sufficient fraction of the universe volume has 
completed the preceding tunneling event. We will show how this condition arises 
through multiple tunnelings, and thereby establish that chain inflation is a viable 
model within string theory or other gravity theories with multiple metastable vacua. 
Although our results are in agreement with the chain inflation models expectations, 
our formulas for the energy density of the vacua, their lifetime and persistence are 
in fact very different from those usually found in the literature (see for instance 
[DISlEllllEllHll?]). In many cases the results have not been derived, but instead 
inferred from calculations based on Old Inflation, not acknowledging that in that case 
only one tunneling step is considered. 

To proceed we write the various contributions to the energy density as 

P^(t) = 5^e„^Pm(t) (9) 

m=0 

m—1 

P^{t) = 5^Ae„^p„(t)^„,„_i(t) (10) 

m=l n=0 

/ m-1 \ 

p^(t) = -AHp' + Y,^^mdtil-Y,Pnit)^m,m-l{t)], (H) 



m=l \ 71=0 



where €m is the value of the potential of the potential energjo in the vacuum m, 
Acm = (cm — Cm-i), Pm{t) IS the fractlou of volume still occupied by the vacuum m 
at time t and, finally, J^m,m-iif) is the energy weighted fraction of uncollided walls at 
time t. The energy density of the ?7i-th vacuum phase is therefore given by 

m 

P^'"* = tmPm + ACm+l ^ Pn(t)^m+l,m(t)- (12) 

n=0 



*Note that in general eo 7^ 0. 



In chain inflation the p„(t) are obtained from a coupled set of differential equations: 

(13) 



Pn 


= -'^nPn 


pN-l 


= —'^n-iPn-i 


Po 


= Tipi, 




>Jp™(t) = l. 



with 



We have 

f^(t) = r^rp^^y«--'(to,t), (14) 

where Tm is the decay rate per unit time and (physical) volume. As we will see later, 
the amount of radiation is negligible and therefore the decay rate will essentially be 
given by the Coleman-Callan tunneling process, and not by thermal nucleation of 
bubbles. 

The uncollided fraction of walls of bubbles of the phase m — 1 is given by (see also 

where Pm,un{t, t') is the fraction of a (reference) wall generated at time t' which is still 
not collided at time t. 

This fraction represents the probability that a point of the (reference) wall has 
not collided at time t, and is equal to the probability that a point just outside the 
wall is still in the phase m, so that: 



Pn 



,{t,t') =e (~'^^'^'" ^« ^'"). (16) 



Note that Pm,un{t', t') = 1, which simply means that at the time of nucleation the wall 
is completely uncollided. 

The solution of the system flT^ is different depending on whether the nucleation 
rates per time Tm are 

(a) equal for every m (that is F^ = T) 

(b) different for every m 

We will show how the correct requirement for having a successful inflationary 
periods can arise, so that chain inflation is a viable model within string theory or 
other gravity theories with multiple metastable vacua. Scenario (b) will be dealt 
with in detail in section 12.41 and we concentrate on case (a) in the following section 

[231 



2.2.1 Percolation and large scale thermalization 

We must now discuss for whicli values of tlie parameters in our model (essentially 
the decay rates) we are able to achieve percolation (and large scale thermalization) 
at every phase transition. There are three time scales that can be used to study the 
final stage of an inflationary expansion. [T2l 113] . 



t™: the time when the volume of the patch starts to contract, that is the time 
when 

phys,m ^^ \ ) 

becomes negative. 

By using Vphys^mif) = o? if)'Pmii) i it is easily seen that f^ is given by the formula 

3i7(a + ^U,j.<0 ^ f „^,^I!^±iM < f„ - 3if (t™) (18) 

Pm Pm\S'c I 

and, in particular Vm — SH^f^) > 0. 

Note that ( fT8l) implies Pm\t=t^ < 0,and therefore t™ > tmax; where t^^x is the 
time for which pm is at its maximum (given by Pm\t=t^^^ = 0). 

t"^: the percolation time for the m-th phase, given by the condition 

total volume of bubbles of phase rn , ^ , , , 

U>OM 19 

total volume '' 

t™ can be found in the following way, [13] . The bubble nucleation rate per unit 
time and coordinate volume is 

.^-^ = a^{t)T^^ie{t-tB), (20) 

where ts represents the time when the nucleation of bubbles starts. 

Since the bubbles expand at the speed of light (in our model) we can relate their 
radius to the time using the formula dr = a{t)~^dt. We choose a reference radius 
ri, corresponding to a time t^i- Bubbles with radius r > ri would nucleate at 
times t < tri- Therefore the nucleation rate per unit coordinate volume for a 
sphere having r > ri is given by 

dT '' 

d^x 



r>n = r^+1 / Ci^t)dt. (21) 



Successful percolation (and, as can be shown, also large scale thermalization, 
see [13]) at time t'T', for ri small, can then be achieved if 

total volume of bubbles with r > ri . Ivr ^dTitZ") ^ „ „ , ,^^. 

\tm > — r'{ — |r>ri ^ 0.34, (22) 

total volume ^ 3 d'^x 

where ^rf is the volume of bubbles of radius ri. 



• t^: the time by which the patch of phase m has become a small fraction of the 
total volume, which for example can be chosen to correspond to the moment 
when p„(Q) = 0.01. 

We will analyze these requirements in more details in the following sections, when 
we will have the explicit form oi Pm{t)- 



2.3 Equal tunneling rate at each step of the ladder 

We will now assume that the tunneling rates are constant and equal at every step of 
the ladder of decay events. That is 

f„(t)~f. (23) 

In this case the solution to (fT3l) is easily seen to be 

p«_i.(«) = -4/- «""■ (2*) 

Furthermore, the uncollided fraction of walls of bubbles of the phase N — k — 1 is 
given by 

where 

r{k + i,t) = k\e-yJ2^^ (26) 

is the incomplete Gamma function. 

It is now possible to determine for which values of the decay rate F successful 
percolation and large scale thermalization are achieved. Formula (1181) for m = N — k 
now reads 

t^-' > ^ (27) 

T-3H, 

and if we take H to be constant during inflatioiij, we see from formula (122!) that 

f logfl - 34^) 
JL(l_e3^(*^i-*™))>0.34 => C> ^^—- ^ (28) 

3x2 3x3, 

where we have used (TT^ to rewrite the formula in terms of F and we have sent tj.^ — ;> 0. 
We can see that if F is sufficiently larger than 3H, then the contraction of volume, 
percolation and large scale thermalization are achieved in very short times. 



^The accuracy of this approximation can be verified by solving the Friedman and Chauduri 
equations using the formula for the energy density of the vacuum that we will provide in (p4|) . 



We can now calculate the vacuum energy, which is given by (see (fT2l) ) 
p^(t) = p^(t)+p^(t) 

(29) 

where we have used the variable y = Tt = -. 

Apart from the last steps of tunneling, it is possible to see that we have 

„V,N-k W,N-k 

^ - ^ -. (30) 



^N-k Aejv- 



k+1 



Note, then, that the transfer of energy from the vacuum to radiation is quite low 
during the main part of the tunneling chain. We can then write: 

N k 

p'it) ~ e-y J2 ^N-k+ij^ + e7ve-^ (31) 



fc=i 



Motivated by the specific model we will describe in section HI we consider two 
possible functional forms for the energy density at every stage 

m^n case I 

m^, „ ^^ , (32) 

^n case 11 

where rrif is a suitable energy scale. In this way we can easily compute the vacuum 
energy with the result 

f;f((A'+l-i)r(iV+l,l) + (i)''«e-i) easel 

l;^(((A'+l-9' + 0r(A'+l.i) + (iV-i + l)(i)~«e-i) oasell 

(33) 
If A^ is large, the incomplete gamma function is to a very high degree of approximation 
constant, T (^N + 1, ^) ~ A^!, for values of t almost all the way up to Nt. We also 

have that (;^) e~~ ^ T (A^ + 1, ^). As a consequence we find, to a high degree of 
accuracy, 

fm^. (A^ + 1-^) easel 
P'^W- J ,. . 12 .. (34) 



^{N+l-^y case II, 



which indeed leads to the expected slow roll. 

It is convenient to define an effective slow-roll parameter 






Using the excellent approximation 

n = N+l--, (36) 

r 

it can be written 

11 I 7; — 77 — case i /__N 

^^-m = { 1 TT • ^37) 

-H^ -^ case II 

For what concerns the number of e-foldings (iVg), a straightforward computation 

gives 

3 3 



iV, = Hdt^ I 



1 n." -n/ y 

3^ 3 ' case I 

71? 

(3^ 



9 9 



2e n 



T- case II 



where we have used H^ ~ ^^p^ and formula (!34l) for the time dependence of the 
energy density. The subscript i/ f indicates quantities at the beginning/end of the 
inflationary period (in particular ni = N using our previous notation). It is easy 
to accommodate the desired value for iVg by choosing nf not too close to rii (but 
not necessarily as low as 1 or 0) and an appropriate value for e. The specific details 
regarding the end of the inflationary period (final metastable vacuum, reheating, etc.) 
will, as we already have stressed, not be treated in this work. 



2.4 Different tunneling rates at each step of the ladder 

In the case when Fj 7^ Tj in (IT5]) . the solution of the coupled set of differential 
equations is given by 

p^_fc = c^- V^--* ( 1 + Yl cf-V(f---f--)A , (39) 

\ i=N-k+l ) 

where the coefficients c^ ~ are given by the recursive equations 

„Af-fc+l _ „N-k„N-k \ / 

''0 — *'0 '^N-k+l ~ 

i N~k+1 

( r^v-fc — ^ N~i ) 

N-k+i ^ ^N-k ^ y ^ f^^^ ^k + 2<i<N (40) 

N-k /f f 

C^-k+l I -*- N-k — i N~k+l 

N 
i=N-k+l 
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As in the previous case we are interested in investigating the conditions for which 
there is a slow-roll regime, and in studying the values for the tunneling rates F^ 
leading to percolation and large scale thermalization. 

In general it is possible to understand how a slow-roll behaviour can arise, by 
observing that every pN-k for fc 7^ is a function peaked (how much depending 
on the Fm's) at a certain time t™*^^. Therefore the tunnelling chain is effectively 
slowed down at every step, until the fraction of volume of the new phase reaches its 
maximum. It is however not possible to write a formula for the total vacuum energy 
such as fl33l) unless a model is specified, so that the precise dependence of the F^ on 
m is known. 

Nevertheless it is possible to estimate the slow-roll parameter as follows. During 
slow-roll we naturally expect: 

H^ ~ ^p^ (41) 



and^' 



From (fT3!) we find 
and from fl32|) we see 



P^ ~ J^ (^mPm (42) 

m 

P^ = - ^ Ae^f „p„ (43) 



Ae^ = — case 1 
m 



(44) 
2e. 



Ae^ ~ case II. 

m 



Then, using this and (HTI) . we find for the slow- roll parameter 



2H 2 ^'"'^^'"^'""' case II 



If we now define an average (-) as 



/_\ = ^m ^m^ mPm 1^ , .„. 

n En ^nPn 



Then the slow-roll parameter is given by 



H f 77Tj{-) case I ,,„, 



^We in general expect this formula because the energy density of uncoUided walls is proportional 
to the energy difference between two consecutive vacua, while the one of the interior of bubbles is 
proportional to the energy level, which is greater than the difference. See also the discussion in the 
previous section. 
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It is interesting (since it will be needed later on, when we will deal with the spectral 
index of perturbations) to compute the time derivative of the slow-roll parameter. 
From fHSl) and using 



e/r\ I 1 fT,„i<^rnT,nPn,m ^ , {J2,n^rnTmPmm l)2j Q^ggJ 



Let us focus on the two terms in each of the brackets. The second one is simply 

(Ene„p„)2 n 

It is the easy to see that using flT^ . the numerator in the first term reads 

Er • -1 — \^r ( —V ^m+i p \ 

C-mJ- mPm ^^ / , J- m+lPm+1 I J- m , ^ -L m+1 I 

^-^ \m m + 1 J 

m m ^ ' 

= -Vf^+iP™+i(Afil^')f„+i + ^A(f„+i)) (50) 
^-^ \ \m + 1 J m / 

where we have defined, for any quantity fm'- 

A iU = frn- fm-1 • (51) 



We find that: 



If we now define 




(52) 



(53) 



Then, from (HHl) and (H7|) . we have 



3i/e2_^^I/^r)(^)-^ easel 

2i7e2_,;2 ^£)-"i_,(£Af )(£)-! casell 



(54) 



where all the averaging has been made using the distribution pm = ^mPm- We note 
that in case I, the derivative of the slow-roll parameter contains a term which will 
also be found in the case of equal tunneling rates at every step (since it does not 
depend on the difference of these decay rates) and a term that is depending on the 
differences themselves. We call them the "universal" term and the "difference" term. 
In case II, instead, we find, beyond the universal and the difference terms, also a 
contribution deriving from the spread of the distribution pm- 
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3 Cosmological perturbations in chain inflation 

3.1 General considerations 

Inflationary models must be able to fulfill a series of requirements such as providing a 
solution to the flatness, entropy, horizon and monopole problems of the standard Big 
Bang. But at the same time, they must also meet experimental constraints coming 
from measurements of the spectrum of cosmological perturbations. 

This last point represents a serious issue for chain inflation, as well as for any 
inflationary model not depending on a single field. Indeed, the CMB analysis strongly 
constrains the amount of entropy perturbations and the form of the spectrum (and 
higher order non-gaussianities) of the adiabatic density perturbations. 

In this section we will derive a set of equations that enable us to discuss the 
cosmological perturbations and compute their adiabatic spectrum. 

3.2 Scalar-type perturbations 
3.2.1 The nature of perturbations 

Within a chain inflation model, density perturbations are generated and evolved 
in a complicated way. Contrary to the usual slow-roll, chaotic inflation model (or any 
model not containing flrst order phase transitions), we have different kind of pertur- 
bations being produced. In a slow-roll model the seeds of cosmological perturbations 
are provided only by the quantum fluctuations that accompany the classical evolu- 
tions of the flelds, metric and fluids. Their dynamics can be described by a set of 
differential equations that relate the perturbations of the different entities (metric, 
flelds, fluids). 

Within a flrst order inflationary scenario, we have instead two kinds of pertur- 
bations being generated. On the one hand the quantum fluctuations as in the case 
of slow-roll, chaotic inflation, and on the other hand the statistical fluctuations of 
bubble nucleation and collision. The presence of these two seeds of cosmological per- 
turbations makes it interesting and important to study and analyze this aspect of the 
model. 

It is clear that we cannot straight away say which of the two mechanism of gener- 
ation is the dominating one. Indeed, the basic quantity we use to describe the energy 
density for vacuum and radiation is represented by the fraction of universe occupied 
by the m— th vacuum 

Pmit) = -r™KP'^y^'^"'(t,t/)p„(t) + r„+iKP'^y^-^'(t,t/)p™+i(t) 

= -Tm(t)pm(t) +Tm+l(t)prn+l{t), (55) 

The quantity Tm{t) represents the average number of bubbles of the m-th vacuum 
generated at a time t, irrespective of the bubbles generated at previous times. The 

13 



actual number of such bubbles is therefore distributed according to a Poisson distri- 
bution with mean value Tmit). Therefore, a first kind of fluctuation in Pm{t) is due 



to the statistical fluctuations in this quantity, being equal to yVmit)- In general, 
solving for pm, we can determine the perturbation in the energy density due to these 
statistical fluctuations. 

But this is not the only seed of fluctuations inpm(t). Also, the volume V^^^'^^'^'^^{t, t') 
depends on the metric which it is measured by. And therefore the quantum fluctu- 
ations in the metric field generate and are sourced by perturbations in the energy 
density of the vacuum and of the radiation as well, due to the coupling between these 
two. 

Our goal will therefore be to obtain a set of equations for the perturbations (relat- 
ing those of the metric and the various energy and pressure densities) , to solve them 
and to obtain the spectrum of the fluctuations. 

We will use the longitudinal gauge, where the most general metric with scalar 
perturbations can be written 

ds"^ = -{l + 24>)dt^ + 2ad^Bdtdx' + ((1 - 2^)5,^ + DijE)dx'dx^), (56) 

with B = E = 0. Note that here Dij = didj - I^.^V^. 

Let us now pause for a moment to discuss the general approach to cosmological 
perturbations. The distinction that is made in cosmology between perturbed and 
unperturbed universes is a fictitious one. The actual universe has a metric of the 
form (l56li 1. and we compare it with a fictitious universe where the metric is a per- 
fect Robertson- Walker one. In our case, pm has a well defined geometrical meaning 
independent of the metric used to compute Tm{t), and has to obey the equation 

Pm = ~^mPm + ^m+lPm+1, (57) 

both in the perturbed and unperturbed case. It is therefore possible to determine 
the perturbations in p^ just by writing F in terms of the perturbed fields and find 
the quantum part (for example in V^^^^^'^''^^(t,t')) and adding the relevant statistical 
perturbation. Similarly, we find the form of jF™"*"^'™ to be determined by the geometry 
also in the presence of perturbations. 

3.2.2 Notation and basic quantities 

As before it is convenient to formally distinguish every phase with an index ttl, and 
associate to each one a vacuum and radiation stress-energy tensor, distinguished by 
£ = V,r as in formulas ([3]). We will also distinguish the quantities in the perturbed 
(real) universe form those in the "background" one by writing the first one with a 
bar, the second one without. Hence we write 

A is the perturbed value 

, (58) 

A is the background 



^We consider only scalar perturbations for the moment. 
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with 

A = A + 5A + 0{{5Af). 

In this formahsm the fundamental equations governing the dynamics of the various 
generahzed fluids are 

V^f^^" = Q'r, (59) 

with, in order for the total stress-energy tensor to be conserved, 

^g^;- = ^ ^Q^;- = ^<5Q:'"^ = 0. (60) 

Let us focus, since it will be specially important in the following, on the vacuum 
sector. First of all, the energy momentum tensor is given by 

Tl. = E('^'''" + p'^n^K + p'^'^'a,., (61) 

m 

where m^ is the proper velocity for the vacuum, and therefore we have 

Second, recall that the formula for the vacuum energy at every phase is given by 

where p^'"^ is the energy in the interior of the bubbles of the m-th phase and p*^'™ is 
the energy stored within the walls between phase m + 1 and phase m. We also know 
that the dynamical evolution of the energy density is 

P — -i-mP + i m+lP - ImP i \^V 



where 7m represents the conversion of the energy of the walls into radiatioi 

By comparing with fl59|) . it is then straightforward to find that the covariant form 
of Ql'"^ is: 

r\V,m -p -Va rriV,Tn , -p -Va rpV,m+l ~ -Va rpW,m 

= Q'^^ul (65) 



where in the last line we have used ( 162|) . Note also that -u^" is the same for all 
components with different m, since it is just the proper velocity of the total vacuum, 
whose energy momentum tensor is (16T!) . 

The perturbations in the sources Q^'™ can be written as [HI [151 [HI fT7] 

,5Q^p'™ = -g^'^^-^Q^.- (66) 

^Q'r = diiQ^'^^v + f'"^), (67) 



^At the microscopical level this would be a series of complicated processes of collisions of walls 
and decays, but we will instead use a well-defined coarse-grained description. 
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where v is the velocity potential, obtained by expanding the spatial components of 
the average proper velocity m* in harmonic functions as in [151 [TB]. In our case, for 
£ = V,r, we find 

^ = E ' .\p v^^- = Y^v-- = v^. (68) 

For ease of notation, we will also use 

5q={j) + P)v, (69) 

where 

5q ^ Y.^q' (70) 

dq^'"^ = (p^'™ + p^.™)^;^''-. (71) 

By expanding ( !65l) in background plus perturbations, writing the proper velocity in 
terms of the velocity potential, comparing with (167]) and using ( l68l) we find for the 
vacuum 

Q^'^'v^ = Q^'^t; + Z^''" = Q^-'^w^ + /^'™. (72) 

In order to study entropy (isocurvature) perturbations, it is convenient to define 
the gauge invariant relative entropy perturbation between fluids i and £', 

Se,e' = -3 (7^, - 7^,0 • (73) 

In the usual definition, for hydrodynamical matter, Si^i' is defined in terms of the 
curvature perturbation on uniform density hypersurfaces. Since we are now studying 
a system where many scalar fields play a role, the correct definition (see [18]) is in 
terms of the comoving curvature perturbationu 

7^, = - Hve. (74) 

The total density perturbation is given by ([T^): 

n = <p-Hv. (75) 

By using these definitions and (168|) . we find: 



f = Y. f"^'^ = ^''(^'' - ^'■) = h^y>r (76) 



3H 



^In |il8j this is indicated as C, which is unfortunately also used for the curvature perturbation on 
uniform energy density hypersurfaces. 
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3.2.3 Equations of motion and spectrum of perturbations 

The most general set of equation for a set of generalized fluids (including scalar 
fields) for the case of pure Einsten-Hilbert theory has been written in [17] (for gen- 
eralized theories of gravity, the enlarged set of equations is in (T6]). In particular we 
will use (see also [14j): 

^ = (77) 

' 3H{H(f) + (f))= ATxG6p (78) 



a 



2 



{H(P + 0) = -47rG(p + p)v = -A-nG5q (79) 

dp'"" + 3i7(5p^'™ + ^P^'") = 3(p^'" + P^'") U + -v^'A + Q^'"^0 + ^Q^''" (80) 

dq'"^ + 3/7(5g^'™ + (p^'"^ + P^'™)0 + (5P^'™ = Q^'™^^ + /^'"^ (81) 

p + p 

(5g + 3i7(5g + (p + P)0 + 5P = (82) 

Note that (5p^'™, (5p^'™, etc, include the contributions from quantum as well as 
statistical fluctuations. In order to re-write this system of equations in a useful way, 
we need to establish the relation between the total pressure and energy densities. We 
already know that 

5P^ = ^Sp^, (83) 

but we also need to flnd the relation between the corresponding quantities for the 
vacuum energy. This can be done by inspecting the relevant equation of motion. 
From 0591) ■ we have: 

Since, according to our discussion above, this relation is valid both for the unperturbed 
and for the perturbed universes, we immediately obtain 

Sp^'"" = SQ^'"^. (85) 

By comparing this result with the general equation ( IHOl) for i = V, we flnd 

f)V,m 
^pV,m ^ _^pV,m ^ ^1^0^ (86) 



3H 



and therefore 






3H 

m 
V,m 



-Sp + Sp^ + J2 ^<P- (87) 

m 
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If we compare this with the result we would obtain in a standard single field slow 
roll or chaotic scenario, we see the importance of the presence of perturbations of the 
radiation at every step of the tunnelling chain. In fact, dp^ plays here a role similar to 
the perturbation of the scalar field kinetic energy in standard slow roll. It is therefore 
important to keep track of it, as it will contribute substantially to the dynamical 
equation of the perturbations. 

On the other hand, from the definition of the momentum perturbation in (171 p . we 
see that 

Sq^^"^ = 0. (88) 

Therefore, using flHTl) for £ = l^, we find 

r)V,m / ^-2 \ 



3H{p + P) V 47rGa2 

and so, summing over all m's: 



sp'' = .JJ, r.A sp+7^<p]+r- (90) 



3H{p + P) V ^T^Ga^ 

Then, from ( IHTl 1901 1761) and neglecting a sub-leading term oc p'', we get 



^"'—E^T^-I^^ + EV^v-,. (91) 



3H ^ AttGu^^ ^^ H 

m m 



where we have used: 



^^H^^QV.m ^ ^^nt_ = 47rGrf 3(2-6) 
3H ^^ 3H 3H dV IQ-kG 



m 



As a consequence 

SP = SP^ + SP' 

A 4 k' ^ I QV,m ^ QV,m 

m m 

Using ( l92l) together with ( 193|) , we finally get the system of equations (from (1791) and 

am) 



(a0)- = 


-AttG a 5q 




6q 


(l ^ ^' 


1 H 


P + P 


y 3 4nGa^{p + P) 

4 H 

3 47rG(p + P) ^'' 


3AnG{p + P) 



(94) 



It is convenient to introduce the new variables TZ, ^ defined by 

a(j) = A-nGHi (95) 

5q 7^ 47rG, , , 

-77 + — e, (96) 



p + p H a 

where we easily recognize in TZ the comoving curvature perturhatiorr^ In terms of 
these new variables, the system becomes 

• _ a{p + p) 

For the moment we concentrate on the adiabatic perturbations, neglecting the 
term proportional to Sy^r- Non-adiabatic isocurvature perturbations (as well as non- 
gaussianitiesLj) are expected to be present in chain inflation models because of the 
presence of various fields that can follow trajectories in field space different from the 
infiaton one, but this depends strongly on the specific features of the various models, 
which our very general picture does not accommodate. 

General arguments justifying the negligibility of the isocurvature perturbations are 
often found in the literature, but they do not apply in our case, or in any other case 
with coexisting fields and conventional hydrodynamical fluids that interact. Lacking 
a general argument regarding the amount of the isocurvature perturbations, we would 
need to rely on specific analysis of the detailed models. This we do not attempt to 
do, and our conclusions in what follows are valid only in those cases were the detailed 
model does indeed contain a negligible contribution of isocurvature perturbations. 

As a word of comment, note that the difference between our final result and 
the usual one (for non-interacting fields and fluids) is, apart from the entropy term, 
proportional to ^^ ~ e (see ( |93|) ). It is therefore, as one could expect, due to the 
interactions. 

In order to normalize the amplitude of quantum fluctuations, we need the action 
governing their dynamics. Our setup, being so general, lacks this information. Nev- 
ertheless, following [12], we can infer the action from the equations of motion (1971) 
themselves, up to a time independent factor. Following the same steps as in [19], we 
introduce the canonical quantization variable <;^ = zR. with 

oxp^rv^i o y 

where we have used the conformal time defined as rj = J a^^dt. The time independent 
operator O can be determined by comparing our result with well-known ones (for 



^"In [T^ this is indicated as C,. 
^^See for example [20] . 
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example with the pure radiation case). We find then that O = — A;^. The action for 
the variable <;^ is in our case found to be (up to a total derivative term) 

S=^l [,'2 + 1,(A + n'e), + '-e] dnd\ (99) 

and the equation of motion is 



with the obvious long wavelength solution <j ~ 2; — i> 7^ ~ const. 

" " \ — I 

Since ^ ~ ^- ~ 2{Ha)'^ we fincO following 



°^ ^^^ Plank ^^^ 

for the spectrum of perturbations, where quantities have to be evaluated at a time 
given by aH = -j=. 



4 A potential for chain inflation from flux com- 
pactiflcation 

In order for chain infiation to be a serious alternative to chaotic infiation, one needs 
a framework where the required potential arises in a natural way. As we will see such 
a framework is provided by fiux compactified string theory. 

4.1 Series of minima from type IIB flux compact iflcat ions 

A particularly interesting compactification of string theory is type IIB with fiuxes. 
It serves as a useful playground for constructing models of the early universe, and 
exhibits in a rather explicit way the problems and promises of the string landscape. 
The potential for a fiux compactified type IIB string is given by 

V {z, t) = e^g^D.iWD^W + e^g^^DrWD^W - 3e^ \Wf , (102) 

where W is the superpotential and K is the Kahler potential. We will focus on the 
dependence on the complex structure moduli and are therefore mostly interested in 
the structure of the superpotential. It is obtained from the fiuxes as 

W = F-U-tH -U, (103) 



^We use the reduced Plank mass A/pi = (SttG) 
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where the period vectors are defined to be 11 = (J^ ^) /g fi, •••)^) with Q as the 
holomorphic three form. 

As shown in |H1 E] it is rather easy to find fiuxes where the potential has a min- 
imum somewhere in moduh space. Of particular interest are minima where one of 
the components of the fiux vector, say F2, is much smaller than the corresponding 
conjugated component, i.e. Fi. We can then make use of conifold monodromies to 
generate new minima. Specifically, we need to move in moduli space around the 
conifold point where the cycle Ai vanishes. When we do that we go through a cut 
and end up on a new sheet. Furthermore, the cycles transform as -Bi — > i?i + nAi, 
where n is the number of times we go around the conifold. The potential after the 
transformation can equivalently be evaluated using the old periods but with 

F(3)^F(3)+n(F2,0,0,...,0) (104) 

i7(3)-^if(3)+n(ff2,0,0,...,0). (105) 

If F2 -C -Fiwe can be reasonably sure that we will find a new minimum near the old 
one on the next sheet. In this way we can generate an infinite series of minima if the 
Fi —>■ 00 limit of F(3) corresponds to a potential with a minimum. The minima will 
be positioned along a spiral staircase turning around the conifold point. 

In [Sj it was discussed in detail how to obtain such series. In [9] it was sug- 
gested that one could make use of geometric transitions between different Calabi-Yau 
manifolds to help generate the series. 



4.2 A potential for chain inflation 

Instead of working with explicit fiux compactifications for specific Calabi-Yau mani- 
folds, we will capture the essential features of the potential with a simple toy model. 
The main results will depend on just a limited number of parameters that in turn 
will depend on the complex geometry of the actual compactification. 

The value of the potential at the minima can effectively be described using the 
results of [S] and general properties of the potential used in type IIB flux compacti- 
fications. There are two cases that are of particular interest. In [S] it was found that 
when the axiodilaton modulus is fixed by the fiuxes themselves the potential at the 
minima has a form such that 

Vn = m'jn = m^ipn Case I. (106) 

If the axiodilaton is fixed by other mechanisms we instead expect 

V„ = -m\n^ = -m^ifil Case II. (107) 

2 -^ 2 

n is in both cases given by the fiux Fi, and we have assumed that we are in a limit 
where the fiux is large but the the minima still remain. In general the potential V is 
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the sum of an underlying m?ip (case I) or ^iv?Lp^ (case II) and superimposed barriers. 
What we have found in this way, is a realization of the chain inflation potential studied 
in the previous section. 

It immediately follows that the energy difference between two adjacent minima is 
given by 

AV^ = m^ case I 

{ (108) 

AV^ = m%n case II 

We can also calculate the field values at the minima with the result 

— rn case 1 

™ ,, (109) 

-^n case 11, 

and the distance in field space between two adjacent minima is given by 

{ml 
-^ case 1 , , 

- (110) 

-^ case 11. 
m 

As we see, there are two energy scales involved, m and mj, that ultimately depend 
on the details of the compactification. irif is given by parameters involving coupling 
strengths and the size of the extra dimensions, while m crucially depends on details 
of the complex geometry determining the distance in field space of the minima from 
a conifold point. In the following we will treat irtf and m as free parameters. 

Let us now turn to the tunneling. In our model the Coleman-Callan decay is 
dominant, hence we have that the decay rate per unit time and physical volume i^'^^' 

-B 



^e"^ = e 2 (^^. (111) 



S= dipW (112) 

is taken over the path that minimizes the integral. In our discussion of chain inflation 
we have assumed a potential depending on a single (real) scalar field. The potentials 
generated from flux compactifications are far more complicated. Even if we restrict 
our attention to just one complex structure moduli, the problem is essentially two 
dimensional with a series of minima winding around the conifold point. We will not 
attempt to perform a complete analysis of the problem, instead we will focus on a 
simple order of magnitude estimate of the ingoing quantities. 

We assume that the least action will essentially be given by a path encircling 
the conifold point with a length given by the circumference. Since, for large n, the 
potential V scales like n or n^, it follows that the barriers between the two minima 
have a height of the order 

I min case I 

I rrirn case 11. 



^•^Recall that the decay rate per unit time is given by r„ — r„ x physical volume. 
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As a consequence we find 



5~Ay,A^/2~H^ "^'"^ (114) 

I -^ n case II, 

and the exponent becomes 

infi^^y^n' easel 

1 277r2 /mf\4 ^^ (115) 

l^(i;f) ^ case II. 

In order for the tunneling not to be too suppressed, and r„ = r„ x physical volume 
to be larger than one, we need that B is of order one. 



5 Matching cosmological data 

We have found that the primordial spectrum is of the forno 

This leads to the spectral index 

1 (A2 )■ e 

""' H Al ^ He 

l-he + jj{lAV){l)-^ easel 

l-4e+i4(£r + i(£Af)(£r casell, ^^^^^ 

n 

where we have used formula (1541) . 

Note that in the case of equal and constant tunneling rate for all tunneling events, 
we get the result 

{1 — 5£: case I 
1 A TT (^^^) 

1 — 4e case 11. 

For case II this is identically the same result as for the case of slow-roll with a quadratic 
potential. The reason for this agreement is that the quadratic potential in slow-roll 
models is special in the sense that it gives rise to dynamics where the scalar field 
rolls with a constant speed. This corresponds perfectly to our case II with constant 
tunneling time r. For other slow-roll models with different potentials, such as the 



""^^This is valid for models where the amount of isocurvature perturbations is negligible compared 
to that of the adiabatic ones. As we have pointed out, there is no general argument to establish 
when this is the case. 
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linear potential of case I, the results based on equal tunneling time will differ from 
ordinary slow-roll based on Hubble friction. 

In general, in order to satisfy the requirements on the spectral index, we expect 

the contributions in (11171) due to a- , and to the difference in tunneling rates to be of 

^ ^ r/n 

order e at most. As a confirmation, it is possible to see that using formula flllll) . we 
would obtain a contribution Be to the spectral index from the dependence on n of the 
exponent. There is also a further dependence on n in the prefactor of the potential 
that is more tricky to derive. We will not attempt to do this. We just note that given 
the slow-roll, the corrections to the spectral index have to be fractions of e. Let us 
use this for a rough estimate of the ingoing parameters. 
Matching the observations require^ 

^l=V, (119) 

where r] ~ 2.5 ■ 10^^. The Hubble constant is then given by 



H= Jsn^ri^Mpi (120) 

which is way below the Planck mass. Writing the Hubble constant in terms of the 
energy density, leads to 

{ (8\/37r^ 77 £:) 4 n~ 4 Mp; easel . , 

^ o 11 1 (121) 

{8y/37i^r]e)-i24n-2Mpi case H ^ ^ 

and therefore ruf is also below the Planck scale. 

To estimate the parameter m we must use the expression for the tunneling ampli- 
tude, (TTT5D . We find that 



2^1 1 



'27- ,12 



716 TUf case I 



m={)i,'M , ' .. (122) 



- 277r 
2B 



■)inimf case II 



Since measurements indicate rig ~ 0.97, we must have £ ~ 10 ^. With e = 10 ^, 
the normalization of the spectrum requires -^ ~ 2.6 ■ 10~^. With n ~ lO'* we find 

h.G-lO-^Mpi easel ,^^^, 

rrif ^ { ^ (123) 

^ 1 9 ■ lO-^Mpi case II 



and with i? ~ 1 we find 



^^We take the data from the 5-years WMAP survey, ^21, 



I 5.3 • lO-^Mpj case I , , 

mr^ { ^ "' 124 

^ 3.1 • IQ-^Mpi case II 
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These values for m and rrif are quite reasonable and can easily be lowered by increasing 
n furtheio. 

It is also of interest to note that the field values at the highest minima are quite 
small, from fll09p 

(2.2-10-^Mj,i easel ,^, 

V^n =~ S a • 125 

[2.7- 10-^ Mpi case II ^ ^ 

Contrary to the case of slow-roll in chaotic inflation, we see that we need not move 
over large regions in field space to achieve inflation. Comparing with chaotic inflation 
one finds that the average motion of the field is in fact much slower than would be the 
case in a potential of the form |m^y9^'^ for the above field values. The reason for the 
slow-roll is the presence of the barriers that make the field temporarily stuck before 
it tunnels. The barriers take the role of the Hubble friction term, which otherwise 
is too weak for slow-roll. To be precise, the average motion of the field if in chain 
inflation corresponds to 0^ ~ e^, rather than a single e as in the slow roll regime of 
chaotic inflation. In chain inflation we instead have a contribution from radiation, 
Pr ~ e, that takes over the role of the kinetic term. We should also add that it is 
straightforward to verify that we obtain the correct number of e-foldings with the 
parameters we are using. 

It is interesting to note that the different speed of sound appearing in the expres- 
sion for A^ effectively increases the strength of the scalar modes compared with the 
tensor modes with a factor three. On the other hand, the chain inflation scenario 
sources additional gravitational waves through the collisions of the bubble walls, be- 
side those due to quantum fluctuations in the metric field (which are the only ones 
present in the usual slow-roll chaotic scenario). These effects constitute an obser- 
vational signature of our model. However, to really verify them, one must analyse 
in details the physics of bubble collisions and derive a precise relation between the 
spectral index and e. We expect this to be extremely model dependent and defer 
such an analysis to future work. 

6 Conclusion 

In this paper we have found that chain inflation can be a viable model of inflatioio. 
We have established that a series of tunneling events can give rise to a slowly changing 
vacuum energy much in the same way as Hubble friction gives rise to slow-roll in 
chaotic inflation. We have calculated the spectrum of adiabatic perturbations and 
found it to be very similar to the one of chaotic inflation up to a difference in the 



-•^^It is possible then to check that with these values, in order to provide the correct number of 
e-foldings and satisfy the requirements of percolation, large scale thermalization and contraction of 
volume at every phase, the inflationary period ends before the field tunnels to the very last minima 
(with vacuum numbers 1 or 0). 

^^In order to establish this in a complete and rigorous way, it is still necessary to study the 
reheating stage and the actual contribution and presence of isocurvature perturbations, which we 
leave for future research, when more detailed models will be available. 
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effective speed of sound. This is due to the fact that radiation plays a similar role in 
chain inflation as the kinetic term for the inflaton does in chaotic inflation. 

We have also found that general properties of the type IIB string landscape can be 
used to generate chain inflation. The key is special properties of flux compactifications 
that naturally give rise to series of minima useful for chain inflation. The relevant 
potentials are not flat or fine-tuned to be flat. Instead, one carefully needs to choose 
the value of other parameters that determine the amplitude for tunneling between 
different vacua. Nevertheless it is interesting to see that one need not have particularly 
extreme values to achieve inflation. 

In order to estimate the likelihood of finding an appropriate potential one needs 
to examine the statistics of the topography of the landscape in a way that has not 
been done so far. We expect this to be a quite challenging task. 
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